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We present a novel type of weighted scale-free network model, in which the weight grows 
independently of the attachment of new nodes. The evolution of this network is thus de- 
termined not only by the preferential attachment of new nodes to existing nodes but also 
by self-growing weight of existing links based on a simple weight-driven rule. This model is 
analytically tractable, so that the various statistical properties, such as the distribution of 
weight, can be derived. Finally, we found that some type of social networks is well described 
by this model. 
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The past decade has witnessed an explosive advance in the understanding of the net- 
work structures emerging in many fields, such as networks of protein-protein interaction, 1 ) 
the WWW, 2 ' and the Internet. 3 ' The most remarkable salient topological feature of these 
networks is scale- freeness, that is, the power-law degree distribution. Theoretical studies have 
revealed two essential mechanisms to generate such scale- free networks: 4 ' growth through the 
continuous addition of new nodes and preferential attachment of new nodes to the existing 
nodes with higher connectivity. Thus, we can say that the scale- free properties of the network 
can be successfully explained by these two mechanisms, though some alternative models with 
the use of quenched disorder fitness distributions have been proposed for generating static 
scale-free networks. 5 ' 

In the above cases, we take into account only the topological network structure, in which 
the links between nodes are either present or not. However, beyond such purely topological 
structures, the interaction strength through the link, i.e., the weight of the link plays a crucial 
role in real-world networks, particularly when we need to consider some dynamical systems 
on the network. In the network of airports, for example, the number of passengers traveling 
between two airports can be regarded as the weight of the link connecting these airports. 6 ' 
Similarly, coauthorship 6, 7 ' and natural language 8 ' are known to be weighted scale- free net- 
works. In these weighted scale-free networks, it is reported that not only the distribution of 
the degree of the nodes but also that of the weights of the links obeys a power law. Hence, 
we will attempt to understand how such weighted scale-free networks appear, in other words, 
whether there exist some specific mechanisms underlying the power law of distribution of link 
weights. 

To account for these power laws, several models have already been proposed recently. 
Most of the previous studies 9-11 ' modeled weighted networks assuming that the weight once 
assigned either remains unaltered or is readjusted only when new nodes are added. Wang and 
Zhang 12 ' reported a model network which grows through preferential attachment. The growth 
is then determined by the fitness and the degree of nodes independently of the weights of links. 
In any case, the weight of a link does not grow by itself independently of the attachment of 
new nodes and degree of the node. It seems that some real-world networks can be explained 
by these models. However, in many other real-world networks, the weight of a link can grow 
spontaneously through a certain weight-driven mechanism. For example, in the coauthorship 
network of the researchers, a node corresponds to a researcher and two nodes are connected 
by a link with a weight. The value of this weight is defined by the number of papers on which 
the two corresponding persons collaborated. In this case, when two persons collaborate again 
in another paper, the weight of the link increases without making new edges. An excellent 
researcher has collaborated with many other researchers on many papers, which implies the 
sum of the weights of the links connecting to the corresponding node (which we call the 
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strength of the node) is very large. In addition, such an important researcher tends to write 
many papers. This means that links connecting to stronger nodes tend to increase their 
weight more rapidly, which is a characteristic of weight-driven preferential attachment. In 
this paper, we propose an analytically tractable model of weighted complex networks which 
grow through the preferential attachment driven by the strengths of the nodes. We show that 
power-law distributions of the degree, weight, and strength can be derived theoretically and 
that these results are confirmed numerically. Moreover, we demonstrate that the networks of 
coauthorship and e-mail can be well explained by this model. 

Before introducing our model, we define some measures to characterize weighted networks. 
First, the connectivity of a network can be expressed by an adj acency matrix aij, whose 
elements take the value 1 if the node i is connected to the node j and otherwise. The degree 
of node i is then defined by ki = Ylf=i a n-> where N is the total number of nodes. In addition, 
the weight of the link between nodes i and j is denoted by Wij . Let us define the strength of 
node i, Sj as Sj = Y2f=i a ij w iji which is the sum of the weight of all the links connecting to 
node i. In this model, we assume that the links are undirected, so that the adjacency matrix 
ciij and the weight matrix Wij are symmetric. 

We present a set of rules for generating the network as follows (Fig. 1). The network 
initially starts with a single node. Rule 1: at each time step, a new node is added to the 
network and a connection is made to one existing node i, where the probability that the 
node i is chosen is proportional to the strength, i.e. Sj/ ^2jLi Sj (strength-driven preferential 
attachment). The weight of this new link is then set to unity. Rule 2: at each time step, ct 
pairs of the existing nodes are selected with the probability proportional to their strengths, 
i.e. SiSj/ (^2,^=1 s fc) • If these two nodes are not connected, they are connected by a link with 
the weight equal to unity. If they are already connected, the weight of the corresponding link 
between them is incremented by one. This rule can be regarded as a generalization of the rule 
in the word web growth. 13 ) Note that the total number of nodes is equal to the time t and 
each node can be labeled by the time u when the node is added. Both the creation of new 
links and the changes in the weight of existing links increase the strength of nodes, and the 
strength of a node increases on average by approximately 2c at each time step. In the case of 
the movie, for example, this implies that an actor/actress plays together with, on average, 2c 
actors/actresses at each time step. 

To analytically obtain the the statistical properties of the network generated by the above 
algorithm, we use a continuous approximation. Now, let us denote the averaged strength of 
the node at time t by s(u, t), where u is the time at which this node was added to the network. 
In the same way as, 13 ) the time evolution of s(u,t) is described by the equation 



ds(u,t) 
dt 



(1 + 2ct) 



s(u,t) 



(1) 



f*dvs(v,t) 
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with boundary condition s(t,t) = 1. Substituting the normalized condition J * dv s(v,t) 
2t + ct 2 into eq. 1, we obtain the solution 



s(u,t) = y/t(2 + ct) 3 1 [u(2 + cu) 3 }. 



(2) 



(2+ct) 3 „- 3 
4 



For <C 1, the distribution of the strength takes the form P(s) 
s(u,t) is approximated by y/t{2 + ct) 3 « _1 2~ 3 . For cu 3> 1, the approximation s(u,t) m 

y/t(2 + ct) 3 u- 4 c~ 3 gives P(s) « | ( t(2 c t^ )3 ) ^ s~ 3/2 ~ l/(2s 3/2 ). Fig. 2 shows the com- 
parison of distribution of strength between theoretical and numerical results, in which the 
exponents obtained in the simulations agree well with theoretical ones. 

Similarly, as a continuous version of the adjacency matrix a^-, let us consider the averaged 
connectivity of the nodes at time t, a(u±,u 2 ,t), where two nodes at each end of the link are 
added at time u\ and u 2 . The connectivity a(ui,u 2 ,t) satisfies the differential equation 



da(ui,U2,t) 

m 



2 ct y(^ [i-a(u 1 ,ti a> t)], 

Uo dvs(v,t)) 



which has a general solution 



a(ui,u 2 ,t) 
1 — exp 



(2+ct) 2 



F(ui,u 2 ), 



\J MlU2(2+CMl) 3 (2+CU2) 3 

where F(u\,u 2 ) is an arbitrary function. Although the boundary condition 

s(u,t) 



a(t, u, t) = a(u, t, t) 



Jo dvs{v,t) 



2+ct 
tu(2+cu) 3 



cannot be satisfied, we set F(u\,u 2 ) = 1 and Taylor expand to obtain 

^ s(u,t) 



a(t, u, t) = a(u, t, t) = 



2+ct 
tu(2+cu) 3 



+ 



So dvs(v,t) ' 

Hence, eq. 3 is an approximate solution of connectivity if F(u±,u 2 ) = 1. 
The average degree of the node born at time u is given by 



k(u, t) 



/ dva(u,v,t) 
Jo 

r2/c r 

/ dv 1 — exp 
Jo 

ft 

1 — exp 



(2+ct) 2 



+ I dv 

'2/c 

r-2/c 



y/uv(2+cu) 3 2 3 
(2+ct) 2 



■^/uv(2+cu) 3 (cv) 3 

= t- J dv exp (-A/y/8v) 



dv exp 



2/c 



-A/ (vW) 



= t 



^ exp (-4*) 



2yfc 



+ f r (0, 4*) 



(3) 
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ex p(-^) + 7¥ erf (^M0 



UttA/^ (V~cA»1) 

where A = (2 + ci) 2 / yfu{2 + cu) 3 and T(a, 6) is the incomplete gamma function. If c is small 
and \fcA <C 1 holds for all u, the degree distribution takes the form P(k) pa yW^^ -3 ^ 2 for 
cu> 1 and -P(fc) pa 9(ct) 3 /(16/c 3 ) for cu <C 1. If c is large, assuming ot > 1 we obtain two 
different degree distributions: P(k) pa ^n/c/k 2 for ^fcA » 1 and P(k) pa y^3/8k^ 3 ^ 2 for 
v^A < 1 (Fig. 3). 

From Eqs. 2 and 4, we find the relationship between the degree and the strength. The 
degree k is proportional to strength s for yfcA <C 1, whereas k pa y^/cs 1 / 2 holds for ^/cA » 1 
(data not shown). The linear relationship for \[cA <C 1 comes from the fact that the weights 
of almost all links between 'young' nodes equal unity. 

As is the case of the adjacency matrix, we can define the continuous version of the weight 
matrix Wij, w(ui,U2,t), whose dynamics are governed by the differential equation 

dw(ui,U2,t) _ s{ui,t)s(u2,t) 

9t ~ (f*dv S (v,t)y 

The solution is given by 

w(u 1 ,u 2 ,t) = (2 + ct) 2 / v / u 1 u 2 (2 + cui) 3 (2 + cu 2 f. 

Note that the relationship J Q * dvw(u,v,t) = s(u,t) is satisfied. The distribution of weights of 
all links in the network is given by 

d(ui,u 2 ) 



P{w) = C f dm 
Jo 



d(ui,w) 



a(u!,u 2 ,t) 



pa C f dui-^-ll — exp(— w)], 
J 2/c dw 

where C is the normalization constant. Using P(w) ~ ^ for large w, we obtain P(w) ~ w~ z l 2 
if cu 2 » 1 and P(w) ~ w~ 3 if cu 2 <C 1 (Fig. 4). In addition, it is often observed that the 
average weight scales with the degrees of the nodes as (u>ij) = (hkj) 6 ' We obtain 9 = 2 for 
the first regime of eq. 4 if cu, L » 1 and cuj » 1, and 9 = 1 if \[cA <C 1 for all u (data not 
shown) . 

The model of the weighted scale-free network we have presented in this paper is simple 
enough to be analytically tractable, which enables us to easily derive the statistical properties. 
In particular, only a single control parameter c determines the network properties, such as the 
distributions of degree, strength and weight. In the coauthorship network of the researchers, 
the quantity c can be estimated as c = 1.5 x 10~ 4 by the condition that, in the real data, t is the 
number of researchers (100945) and 2t + ct 2 must be equal to the summation of the strength, 



5/?? 



J. Phys. Soc. Jpn. Letter 

1.75 x 10 6 . The smallness of c implies that new papers on which no new researcher collaborates, 
is rare. The network generated by this model with the above estimated c exhibits scale-free 
properties similar to the real coauthorship network of researchers (Fig. 5). The important 
point is that the various scale-free properties and the exponents stem from a single real- 
measured parameter c. The actor /actress collaboration network can be fitted quite well using 
the present model (data not shown), and this network also has small c. On the other hand, 
the e-mail network can be regarded as a typical case of the large c, because an enormous 
number of e-mails are communicated everyday, regardless of whether new persons begin to 
use e-mail or not. This means that the weight of the link between the existing nodes tends 
to increase independently of the addition of new nodes. The exponent of P(k) the of e-mail 
network is reported to be around 1.8, 14 ** which lies between the exponents 3/2 and 2 of our 
model network with large c. 

In conclusion, we have proposed a novel type of weighted scale-free network model. The 
significant, novel rule in this model is that the weight of the links associated with some 
constant fraction of existing nodes (represented by the parameter c) spontaneously increases 
independently of the attachment of new nodes. As a result, two types of scale-free network 
emerge depending on the parameter c. The resultant networks for small c and large c seem to 
capture the statistical properties of the coauthorship network and e-mail network, respectively. 
This suggests that the proposed simple algorithm is suitable for studying certain types of real- 
world social networks. Another important point is its analytical tractability, which means that 
some statistical properties can be derived theoretically in this model. This is helpful not only 
for a deeper understanding of the weighted scale-free networks, but also for developing some 
extended version of the model, and thus for studies in the other fields related to complex 
systems, such as oscillators, 15 ) epidemics, 16 ) and biological networks. ^ We believe that this 
model provides some insights on the dynamical evolution of such social networks, leading to 
the understanding of more general mechanisms underlying complex networks. 
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Fig. 1 



At each time step, a new single node (a blue circle) appears and connects to one existing node 
with a link of weight one (a blue link). This new link is created by preferential attachment 
with the probability proportional to the strength of the existing node. At the same time, some 
pairs of existing nodes are chosen on a simple strength preferential rule (see the main text 
for details), and the weights of the links between these chosen nodes (a green link) increase 
by one. If no corresponding link exists, a new link of weight one (a red link) is created. The 
numbers on the nodes and near the links indicate the strengths and the weights, respectively. 



Comparison of distribution of strength between theoretical and numerical results. Note that 
for the network with c = 1.75 x 10~ 4 the power law exponent changes at the crossover point 
indicated by the arrow. The bin width is set to 1 for s < 100 and 100 for s > 100 because 
points are sparse in the region s > 100. 



work c = 1.75 x 10~ 4 , P(k) « ^/3/8k~ 3 / 2 (k < k c ) and 9(ct) 3 /(16/c 3 ) (jfe > k c ) (k c = Sir /(3c)). 



Distribution of weights. No crossover behavior is observed, because cu <C 1 (cu S> 1) holds 
for almost all nodes in the network c = 1.75 x 10 -4 (c = 0.5). 



Comparison of scale-free properties between the coauthorship network (filled circles) and 
the present model (cross). Degree distribution (top left), strength distribution (top right), 
weight distribution (bottom left), and strength-degree relationship (bottom right) are shown. 
The coauthorship network is reconstructed from the Geological Literature Search System 
(GEOLIS+ CD-ROM Ver.5) provided by AIST (permission number 63500-A-20070322-001). 



Fig. 2 



Fig. 3 




Each point of crossover k c is indicated by an arrow. 
Fig. 4 



Fig. 5 
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